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Abstract 

In the framework of perturbative QCD in the fourth order over the ag couphng 
constant for the production of two pairs of quarks and in the model of weak 
binding of the quarks into meson, the analysis of the hadronic production of 
mesons, containing 6-quark, is performed, and a minimal transverse momentum 
is determined, so that at the momenta greater than the found one, the meson 
differential spectrum can be reliably described in the model of factorization of 
the hard 6-quark production and the subsequent fragmentation into the mesons. 
At low transverse momenta of the meson, nonfragmentational contributions are 
essential. They are determined by the complete set of QCD diagrams in the given 
order over q^, so that the latters result in the effect of destructive interference 
with the diagrams of fragmentation at z = 2E/^/s close to 1. 



1 Introduction 

In the study of the heavy quark production mechanisms, the consideration within 
the perturbative theory is basic. This is related with the fact that the heavy quark 
mass rriQ determines the scale being much greater than the confinement scale in QCD, 
rriQ ^ ^QCD- Then the cross-section of the heavy quark production is factorized as 
the partonic cross-section of the hard production of QQ quarks and the distribution 
functions of quarks and gluons, in the interaction of which the QQ pair is produced. 

Next step is connected with the heavy quark transform into hadrons containing 
these quarks. Ordinary in this consideration, one uses the model of the heavy quark 
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fragmentation into the quarkonium Q {Qq) + X with the fragmentation functions 
measured in the (Qg)-hadron production in +e^-annihilation, for instance. Such ap- 
proach is the generahzation of the factorization theorem on the case of hadronic pro- 
duction Apparently, the scheme mentioned above does not satisfactory work as one 
can see, say, in the 5- meson production at the FNAL colhder energy p|. Indeed, the 
result of theoretical predictions for the 6-quark production in the 0(ag)-calculations 
is in the systematic disagreement with the experiment. On the other hand, the pro- 
duction of mixed flavor heavy quarkonium. Be, considered in refs. p, ^ ^ , where the 
hadronization is already included as the composite element of the model, points out the 
deviation from the naive intuitive expectation of the fragmentation regime dominance. 

The detailed consideration of the i?c-meson production processes shows, for in- 
stance, that in a broad kinematical region these processes can not be reduced to the 
simple consequence: 1) the 6-quark production, 2) the 6-quark fragmentation into the 
i?c-mesons with the fragmentation function D^^bX'^) known from the calculations in 
e"''e~-annihilation. So, for the B^. production in 77-collisions considered in 0, the 
c-quark fragmentation into Be is greater than the b ^ Be fragmentation. Remember, 
that in e"'"e~-annihilation the c Be fragmentation is suppressed by two orders of 
magnitude with respect to the b ^ Be fragmentation [0. 

In the hadronic production, where there is no such enforcement of the c-quark 
contribution due to the [Qe/QbY = 16 factor, the 6-quark fragmentation into Be is 
still not dominant. The naive expectations of the fragmentation dominance are just 
reahzed at px > Mb^ (pt > 30 GeV for Be and pr > 40 GeV for B*). By the way, the 
contribution of diagrams different from those of the fragmentation type on the topology 
(the recombination contribution in the rest of the paper) is basic. 

In this paper we investigate in what measure these additional contributions can be 
essential in the 5-meson production. We evaluate the transverse momentum, so that at 
momenta greater than the determined one, the heavy meson production can be reliably 
described by the fragmentation model. We study these problems in the framework of 
model for the production of weakly bound quarkonium with the following standard 
assumptions. 1) The calculation of the production cross-section for two pairs of quarks 
is performed in the leading Born approximation (diagrams of the fourth order over a^). 
2) We neglect the binding energy and relative quark motion inside the meson. The 
latter means that the quarks inside the (Qg)-meson move with the parallel momenta 
Pq = Pq~- Since the derivation of analytic expressions for the considered processes 
looks to be problematic, we use a numerical way for the evaluation of amplitudes 
corresponding to 36 Feynman diagrams in the hadronic production and to 20 diagrams 
in the photonic production of quarkonium. 
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2 Factorization of inclusive spectrum of heavy me- 
son 

The basic subject of our discussion below will be the factorization of the inclusive cross- 
section for the hadronic production of heavy meson in the region of high transverse 
momenta. In present Section we consider the subprocess of the gluon-gluon production 
of the meson. The physical basis for the construction of the factorization model is the 
appearance of two energetic scales determining the hard production of quarks at high 
transverse momenta and a soft nonperturbative forming of the bound state. The case, 
when the conditions i/i ^ nig ^ nig ^ Aqcd are valid, is of a special interest. 
At the given ordering of scales, the quark production allows the description in the 
framework of perturbative QCD, and their hadronization into the meson is described by 
the nonrelativistic model of weakly bound state: the energy of the quark binding is low, 
^Qq ^q, nT'Q, ^ud the perturbative production of quarks takes place in the region, 
where the quarks composing the meson move with the same velocity equal to the 
meson velocity. Then the quark momenta are equal topq = uiq-v^ , = irtq-v^. The 
mg 3> mq condition provides the appearance of two scales in the production process 
of quarks with different flavors, so that we can talk on the Q-quark fragmentation into 
the (Qg)-meson. 

The technique of such calculations was given in refs. [0, |, ^, in details. 

As in the model of fragmentation, in the 0(ag)-order over the QCD coupling con- 
stant the cross-section of the subprocess gg {Qq) + Q + q is proportional to the 
ag|\E'(0)p value, which determines an overall normalization of the cross-section and it 
includes all numeric effects of large distances (Aqcd and the meson size < rgq >). In 
concrete calculations we have supposed mQ/mq ~ 17, y/s/niQ ^ 20, niq/ Aqqb — 10. 

By the general theorem on factorization [|T| it is clear that at high transverse mo- 
menta the fragmentation of the heavier quark Q — > (Qq) + q must dominate. It is 
described by the factorized formula 



where /x is the factorization scale, da/dhx is the cross-section for the gluon-gluon pro- 
duction of quarks Q + Q, D is the fragmentation function. In what follows, we consider 
the production of the S'-wave states in the (Qq) system: the vector 1^ state and pseu- 
doscalar 0~ one. For these states the scaling analytic expressions of the fragmentation 
functions were derived at M^/s <^ 1 in the leading order of perturbative QCD 0, H, |^. 

The distributions over the meson transverse momentum in the gluon-gluon produc- 
tion are shown in Fig. |I], in which one can see that at > ^ ^{Qq) exact 
perturbative calculations of the 0(Q;^)-order prove the validity of the factorization the- 
orem, and they fix the scale, from which the production mechanism reaches the regime 
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of fragmentation. Note, that numerically at the chosen ratios of the quark masses and 
the total energy of the partonic subprocess we have p™"^ ~ (5 6) • uiq, and the latter 
value is greater for the vector state. The analysis shows also that the low boundary for 
the applicability of the factorization approach shifts to the region of lower momenta at 
the decrease of the mq/niQ ratio. 

Compare now the result of the exact perturbative 0(ag)-calculation for the dif- 
ferential cross-section da/dz, z = 2E/\/I, with the result of the factorization model 
analogous to (0). Remember, that in e"'"e~-annihilation the distribution under con- 
sideration coincides with the fragmentation function. One can see in Fig. |^, that, in 
contrast to the fragmentation model, there is a sizable additional contribution of the 
recombination type in the region of low z (i.e. at low energies and, hence, low pt)- We 
have found that the recombination contribution is generally given by the diagrams of 
the evolution type for the gluon splitting into the qq or gg pairs. However, the advan- 
tage of the given perturbative approach at ^ Aqcd is the exact, process dependent 
account for the finite values of masses, virtualities and transverse momenta, in contrast 
to the universal approximate method of Altarelli-Parisi. 

We have in details considered the contributions of each diagram in the region of 
z — > 1. In the covariant Feynman gauge the diagrams of the gluon-gluon production of 
Q + Q with the subsequent Q {Qq) fragmentation dominate as well as the diagrams, 
when the qq pair is produced in the region of the initial gluon splitting, are large. 
However, as one can see in Fig. |^, the contribution of the latters diagrams leads to 
the destructive interference with the fragmentation amplitude, and this results in the 
"reduction" of the production cross-section in the region of z close to 1. So, /S.a/a is 
about 50%. In the axial gauge with the vector = Pq this effect of the interference has 
still a more bright appearance, since the diagrams like the splitting of gluons dominate 
by several orders of magnitude over the fragmentation, but the destructive interference 
results in the cancellation of such extremely large contributions. By the way, we have to 
note that discussions and speculations given in ref. are certainly misleading, since 
the authors draw conclusions by the consideration of the singularities in the virtual 
particle propagators, which could determine a dominating contribution. As we have 
just seen, some separately large contributions can destructively interfere to cancel each 
to other. This interference is caused by the nonabelian nature of QCD, i.e. by the 
presence of the gluon selfaction vertices. The comparison with the abelian case will be 
in details given in next Section. 



3 Gluon production of Be 

The model for the production of the mixed flavor heavy quarkonium, as it stands in 
previous Section, is the most realistic in the description of the differential spectra for 
the mesons |TD|, for which, however, the rric/mb ratio can not be yet considered as 
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the small parameter. 

As it was shown in ^ ^ and one can see in Fig. ^, the fragmentation regime 
for the 5c- mesons is strongly delayed into the region of high transverse momenta, 
pmm ^ 35^40 GeV. The last fact is related with the comparatively large value of the c- 
quark mass. Thus, one can again talk on the confirmation of the factorization theorem 
for the hard gluonic production of 6-quarks at high transverse momenta pt 3> M^^ 
and the less hard fragmentation of 6-quark into i?^*)-meson {rric -C pr)- However, the 
large numeric value of p^'" points out the fact that the basic amount of events of the 
hadronic i?^*^-production do not certainly allows the description in the framework of 
the fragmentation model. This conclusion looks more evident, if one considers the 
Sc-meson spectrum over the energy (see Fig. |). 

The basic part of events for the gluon-gluon production of Be is accumulated in the 
region of low z close to 0, where the recombination being essentially greater than the 
fragmentation, dominates. This is the significant difference between the Be production 
and the (Q^)-meson one, considered in previous Section, where both the fragmentation 
and recombination have given comparable contributions. As it stands for mg 3> m^, 
one can draw the conclusion on the essential destructive interference in the region of z 
close to 1, for the pseudoscalar state. 

To stress the role of the interference diagrams related to the nonabelian selfaction 
of gluons, we have considered the process with abelian currents. To avoid effects of 
an enforcement for the fragmentation contributions by different quarks due to the 
difference between their charges 0, we have supposed the abelian charges of produced 
quarks to be equal each to other. As one can see in Fig. where the spectrum over 
the energy of the pseudoscalar state is shown, in the abelian case the effect of the 
destructive interference due to the additional contribution by the selfaction of gauge 
quanta, is absent. So, the agreement between the factorized model of fragmentation 
and the exact perturbative calculation is quite good at z close to 1. 

The direct verification of the given mechanism for the i?c-meson production could 
be the comparison of the i?c-meson spectra in two semispheres in the region of the gluon 
fragmentation and in the photon one, in the photonic production of B^ at nucleons. 

4 Hadronic production of B-mesons 

The model of the gluon-gluon production of quarkonium, as it is described in Section 
2, can be, after some notes, expanded to the hadronic production of mesons with 
a heavy quark Q {mq ^ Aqcd) and a light antiquark q {nig ~ Aqcd)- Then one 
can talk on the model of meson composed of the constituent quarks, whose motions 
inside the meson are neglected, and, hence, one can apply the nonrelativistic model of 
quarkonium with weakly bound constituents. As for the nonperturbative parameters 
of model, rUg and |\I'(0)|, they can be phenomenologically fixed over the spectrum form 
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and total cross-section of the 5-meson production in e^e~-annihilation at the Z-boson 
peak, where the fragmentation regime dominates certainly. Considering the a~^da/dz 
spectrum over z = lEj^^ one has to take into account the perturbative evolution 
breaking the scaling for the function of the 6-quark fragmentation into the 5-meson, 
D'^^^{fi = ^/s■, z). After that, the niq = 0.3 GeV mass value agrees with the observed 
form of the fragmentation function, still the accuracy of such estimate is low because 
of large experimental errors in parameters. So, in the one-loop approximation for the 
perturbative evolution the average fraction of energy taken by the 5-meson, equals 

\\n{y/s/AQCB)J 

where n = 32/(96), b = 11 — 2nj/3, nj = 5, and {z{mh)) is the scaling value of the 
momentum fraction determined over the b B^*^ fragmentation functions derived in 
[0, |, 0. From (|^), the experimental value 



{zBimz)) = 0.708 ± 0.003 ± 0.015 



and Aqcd = 85±20 MeV corresponding to ^^(m^) = 0.117±0.005 [jril in the one-loop 




approximation, one finds rUg = 0.2 4-0.3 GeV. In the following calculations we suppose 
rrig = 0.3 GeV. As for the normalization of the total cross-section of the 5-meson 
production at the Z-boson peak, it is determined by the a^|\&(0)p value, where ^1^(0) is 
expressed through the leptonic constant of the 5-meson in the limit of static 6-quark, 

f Stat 
J B 5 

1^ 

In various papers [|12| one gives estimates corresponding to f^'^^ = 200 4- 320 MeV. 
Supposing the relative yields of B^, B^ and i?°-mesons to be equal to 0.4, 0.4 and 
0.2, we get = 0.8 ± 0.2. Of course, this phenomenological value of QCD coupling 
constant is large for one can talk on the reliability of perturbative calculations of the 
i?-meson spectrum in the model of Section 2. Nevertheless, one can believe that, being 
phenomenological, such model is able to describe the basic differential characteristics 
in the hadronic production of B-mesons, as it is in e"'"e~-annihilation. 

In Section 2 we have chosen the quark mass ratios such that they correspond to 
niQ = nib = 5.0 GeV, nig = 0.3 GeV and = 100 GeV. The only difference between 
the real case of the i?-meson production and the consideration given in Section 2, is 
the change of the overall normalization of the total cross-section. The ratios of spectra 
over the B* and 5-meson transverse momenta are presented in Fig. ^, where they 
are calculated in the model of fragmentation and in the 0{aj)-OTdeT of QCD in pp- 
collisions at i/i =1.8 TeV with the use of gluon distributions given in ||T^. We see 
that the cross-section has the significant additional contribution of recombination in 
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the region of px < 15 GeV, that quahtatively agrees with the experimental data from 
CDF [| (see Fig. |). 

One has to note that the gluonic luminosity of the gg B + b + q process en- 
forces the contribution of the region near the threshold, ^/sth = 2(mQ + rriq), where, 
strictly speaking, the fragmentation approach for the 6-quarks at low momenta has un- 
certainties of the principal character, since the physical phase space of the subprocess 
is suppressed in comparison with the two-particle phase space of free 6-quark pair. The 
enforcement of the near-threshold region due to the luminosities and phase space leads 
to the fact that the fragmentation model evidently gives overestimating values for the 
differential cross-sections at low px- 

The uncertainty with the phase space decreases rapidly, if one limits the total energy 
of the gluon subprocess from the down. So, at \/§ > 60 GeV one can clear see in Fig. 
1^, that the fragmentation picture of the 5-meson production does not correspond 
to the exact perturbative 0(a^)-calculations of the meson spectrum in the model of 
hadronization for the weakly bound constituent quarks, at pt < 25 GeV. It is clear that 
at high energies, the multiple production of light quark pairs can take place in addition 
to the single one. The large value of the phenomenological constant as ~ 0.8 points 
out the latter possibility. In the fourth order over as this means that the invariant 
mass of the (bq) pair can be large in accordance with the (qq) pair mass can be close 
to the threshold. Such contribution could make softer the 6-quark spectrum, i.e. it 
could result in the additional contribution into the region of low pt < 15 GeV, where 
the fragmentation regime does not still work. 

Thus, we have shown that the perturbative 0(a^)-calculations of QCD, accord- 
ing to the reasonable model of the 6-quark hadronization, confirm the applicability 
of the factorization theorem at p^ > 15 GeV. They point out the appearance of ad- 
ditional, nonfragmentational contributions at low transverse momenta of the meson, 
p^ < 15 GeV, so that these corrections appropriately reflect the observed experimental 
disagreement with the fragmentation model. 

5 Conclusion 

In the framework of perturbative QCD up to 0(af)-terms and in the simplest model 
of quark hadronization, we have analyzed the differential cross-sections of the heavy 
mesons over the transverse momenta and the energies. This analysis shows the follow- 
ing. 

1) As one could expect by the theorem of factorization for hard processes, in the 
region of Pt > p™™ ^ ^{Qq) fragmentation contribution dominates. At px < p™"^ 
the quark recombination into the heavy meson is basic, so that its relative value depends 
on the mq/mg ratio. One has p™™ ~ 40 GeV for the i?c-niesons, and the recombination 
contribution in the total cross-section is much greater than the fragmentation. One 
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has p™™ ~ 15 GeV for the B-meson production in hadron colhsions, and the relative 
contributions by both the recombination and fragmentation are comparable with each 
other. 

2) The exact perturbative 0(a^)-calculations point out the nonfragmentational con- 
tributions, which were not yet taken into account elsewhere previously, and they dom- 
inate in the region of low z = ^Ej^fs. These contributions increase by their relative 
magnitude in the total cross-section versus the decrease of the mq/mq ratio. 

3) In the region of low transverse momenta and at z close to 1, there is the destruc- 
tive interference of the fragmentation type diagrams with the diagrams of the initial 
gluon splitting into the qq pair, so that the presence of the interference effect is deter- 
mined by the nonabelian selfaction of gluons in QCD, and the effect is absent for the 
abelian currents. 

This work is in part supported by the grant of Russian Foundation of Fundamental 
Researches, A^^ 96-02-18216. 
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Figure Captions 



1. The distributions over the transverse momenta of pseudoscalar and vector (Qq)- 
states, described in the 0(a;f )-order of QCD (the dashed and soUd hnc his- 
tograms, correspondingly) and in the fragmentation model (dashed and solid 
curves) in gluon-gluon collisions. The cross-section is given in arbitrary units, 
the transverse momentum is in masses of uiq/B. 

2. da/dz for the {Qq)-system (notations as in Fig. 1). 

3. da/dpT for the S^*)-mesons at = 200 GeV (notations as in Fig. 1). 

4. da/dz for the i?c-niesons at y/S = 100 GeV. The dashed histogram presents the 
gg ^ Bc + b + c process, the dotted one is the abehan case. The curve shows the 
result of the fragmentation model. The cross-sections are normalized over the 
cross-section of the bb pair production. 

5. The ratio of spectra for the summed contributions oi B* and S-mesons in pp- 
coUisions at ^/i = 1.8 TeV, as they are calculated in the 0(Q;f)-model of the 
weakly bound quarkonium and in the fragmentation model, a / a^'^^fJpx) is 

given as the solid histogram, the dashed histogram is the same but with the cut 
off > 60 GeV. 

6. The CDF data on the differential cross-section of the S-meson production at the 
FNAL Tevatron, in comparison with the 0(Q;^)-prediction of perturbative QCD 
for the 6-quark production, cXexp/crQCD- 
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